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We investigate the finite temperature properties of an ultracold atomic Fermi gas with spin popu- 
lation imbalance in a highly elongated harmonic trap. Previous studies at zero temperature showed 
that the gas stays in an exotic spatially inhomogeneous Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) 
superfluid state at the trap center; while moving to the edge, the system changes into either a non- 
polarized Bardeen-Cooper-Schriffer superfluid {P < Pc) or a fully polarized normal gas (P > Pc), 
depending on the smallness of the spin polarization P, relative to a critical value Pc. In this work, 
we show how these two phase-separation phases evolve with increasing temperature, and thereby 
construct a flnite temperature phase diagram. For typical interactions, we flnd that the exotic FFLO 
phase survives below one-tenth of Fermi degeneracy temperature, which seems to be accessible in the 
current experiment. The density proflle, equation of state, and speciflc heat of the polarized system 
have been calculated and discussed in detail. Our results are useful for the on-going experiment at 
Rice University on the search for FFLO states in quasi-one-dimensional polarized Fermi gases. 
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I. INTRODUCTION 

Impressive experimental progress has occurred recently 
in the field of ultracold Fermi gases. In particular, ex- 
perimental realizations of the crossover from a Bardeen- 
Cooper-Schrieffer (BCS) superfluid to Bose-Einstein con- 
densate (BEC) using a well-controlled Feshbach reso- 
nance [J, 13, BH, S, m 0] have opened an intriguing op- 
portunity to study long-standing many-body problems in 
condensed matter physics, such as high temperature su- 
perconductivity. Most recently, two experimental groups 
at Rice University and MIT have successfully manipu- 
lated a two-component atomic Fermi gas of lithium atoms 
with unequal spin populations [1, 0, [HI, • This 

type of matter is of great interest, and stimulates intense 
efforts on studying an unsolved problem in condensed 
matter and particle physics: what is the ground state of 
a spin-polarized Fermi gas with attractive interactions? 

Without spin polarization, the answer is known, given 
flfty years ago by Bardeen, Cooper, and Schrieffer. Due 
to attractions, atoms of different spin species at the same 
Fermi surface with opposite momentum form Cooper 
pairs characteristic of bosons, and thus undergo a BEC- 
like superfluid phase transition at a sufliciently low tem- 
perature. This BCS mean-fleld picture is very robust and 
is valid not only quantitatively at weak-coupling, but also 
qualitatively in the strongly interacting limit [3], where 
the fluctuations of Cooper pairs become important. The 
key ingredient of the BCS pairing is the fully overlapped 
Fermi surfaces, which create a maximum for the phase 
space. In the presence of spin polarization, however, the 
two Fermi surfaces are no longer aligned. The standard 
Cooper pairing scheme is thus not applicable. For a small 
number of unpaired atoms, or, a small spin polarization, 
some non-standard pairing scenarios have to be devel- 



oped. For a large spin polarization above threshold, these 
unpaired atoms will eventually destroy the coherence of 
pairs and hence the superfluidity. 

The determination of non-standard pairing scenarios 
and the related exotic superfluidity Hes at the heart 
of polarized Fermi gases. A number of pairing pro- 
posals have already been s ugg ested in the literature, 
including breached pairing Il4ll or Sarma superfluidity 
[H, 01, phase separation [l7|, deformed Fermi sur- 
face and spatially modulated Fulde-Ferrell-Larkin- 
Ovchinnikov (FFLO) states [H, [l^l. At the moment, 
a theoretical consensus on the true ground state of a 
three-dimensional spin-polarized Fermi system is yet to 
be reached [2l|, HI, 111, HI, [H, H^l • Experimental ob- 
servations near a Feshbach resonance are not so helpful 
as one would expect, due to the presence of the harmonic 
trap. To take it into account, the local density approxi- 
mation has been widely used [H, [H, \M, Ell, S [H, IH , 
as well as the mean-fleld Bogoliubov-de Gennes equation 

Among various pairing schemes, the study of FFLO 
phases has the longest history — more than four decades 
li^l- In this scenario, the two mismatched Fermi 
surfaces have been shifted by an amount in momentum 
space [19], in order to have a small overlap of the surface. 
As a result. Cooper pairs with a finite center-of-mass 
momentum may form, and thereby support a spatially 
inhomogeneous superfluidity [23|. However, due to the 
much reduced phase space for pairing, the window for the 
appearance of FFLO states turns out to be small in three 
dimensions [2l|,l23|- This makes the experimental search 
for the FFLO states extremely challenging. Till now only 
indirect evidence has been observed in the heavy fermion 
superconductor CeCoIns (s^l- No FFLO signal has been 
found in ultracold polarized atomic Fermi gases. 



2 



Luckily, the FFLO states are theoretically found to be 
favorable in low dimensions, which can be realized ex- 
perimentally using an optical lattice In one dimen- 
sion (ID), where the whole Fermi surface shrinks to two 
Fermi points, the reduction of pairing phase space is less 
significant. Consequently, at zero temperature the FFLO 
phase becomes much more robust O Ell, [ill, compared 
to its 3D counterpart 1221 . This expectation was indeed 
found in recent works |43l. 0, HI, [H, 113, Further, 
the presence of a harmonic trap will not change the essen- 
tial picturOj as firstly shown by Orso [i^ and the present 
authors [46, 47]]. Nonetheless, the trap leads to phase 
separation [45 . |46j : While the system remains always in 
the FFLO phase at the trap center, at the trap edge, it 
can be either a fully non-polarized BCS state or a fully 
polarized normal state, for a small or large spin polariza- 
tion, respectively. Most recently, the spatial distribution 
of pair correlation functions and the momentum distri- 
bution of a spin-polarized ID Fermi gas in lattices have 
been investigated by using numerically accurate density- 
matrix renormalization-group methods [i^, [H3, HH and 
quantum Monte Carlo simulations [H^ . 

Inspired by theoretical suggestions, an experiment 
with a polarized Fermi gases in a highly elongated har- 
monic trap is now underway at Rice University [H3|. 
However, the experiment will necessarily be carried out 
at a finite temperature. It is then desirable to under- 
stand how the two phase separation phases found at zero 
temperature evolve as temperature increases, and in par- 
ticular, what is the temperature window for the presence 
of FFLO states, compared to the lowest experimentally 
attainable temperature. 

In a previous study we have presented a system- 
atic study of zero-temperature quantum phases in a one- 
dimensional spin-polarized Fermi gas. Three theoreti- 
cal methods have been comparatively used: the mean- 
field theory with either an order parameter in a single- 
plane-wave form or a self-consistently determined or- 
der parameter usinp the Bogoliubov-de Gennes equations 
IHtI . [Hi , as well as the exact Bethe 
We have found that a spatially in- 
homogeneous FFLO, which lies between the fully paired 
BCS state and the fully polarized normal state, domi- 
nates most of the phase diagram of a uniform gas. The 
phase transition from the BCS state to the FFLO phase is 
of second order. We have also investigated the effect of a 
harmonic trapping potential on the phase diagram, and 
find that in this case the trap generally leads to phase 
separation as mentioned above. We finally investigate 
the local fermionic density of states of the FFLO phase. 
A two-energy-gap structure is shown up, which may be 
used as an experimental probe of the FFLO states. 

In this work, we address the urgent finite-temperature 
problem by extending our previous zero-temperature 
analysis. In particular, we focus on the use of weak- 
coupling Bogoliubov-de Gennes (BdG) theory [H, [13, 
H|. Based on this mean-field approach we are able to 
obtain the density profiles of the system at finite tem- 
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Figure 1: (color online) Finite temperature phase diagram of 
a ID trapped spin-polarized Fermi gas at a given coupling 
constant 70 = 1.6, as a function of the spin polarization, 
P = (W| - iVi)/(iVt + iV|), where iVt and iV^ are the number 
of spin-up and spin-down atoms, respectively. TV is the Fermi 
temperature of an ideal, non-interacting Fermi gas in a har- 
monic trap. The solid line gives the boundary of the second 
order phase transition from a superfluid to a normal state. 
The dashed and dot-dashed lines distinguish the different su- 
perfluid phases: a pure BCS phase, a phase separation with a 
FFLO at center and a BCS outside (FFLO-BCS), and a phase 
separation with FFLO and normal components (FFLO-N). 
Due to the presence of harmonic traps, these boundaries may 
be better understood as crossover, rather than phase transi- 
tion lines. 



peratures, as well as its thermodynamic properties, in- 
cluding the entropy, energy and specific heat. 

Our main result, a finite temperature phase diagram, 
is shown in Fig. 1 for a typical interaction strength. At a 
finite but low temperature, in addition to the two phase 
separation states mentioned earlier, two new phases — a 
pure BCS state and a partially polarized normal state — 
enter the phase diagram, respectively, at lower and higher 
spin polarizations. The space for the phase separation 
states shrinks with increasing temperature and vanishes 
completely at one-fifth of the Fermi temperature Tp. 
Therefore, even at a temperature as high as O.lTp, it 
is still possible to observe the long-sought FFLO state 
at the center of the harmonic trap. This suggests that 
the FFLO state in ID polarized Fermi gases is indeed 
attainable with current experimental techniques 

In the following section, we briefiy review the mean- 
field BdG theory of the ID spin-polarized Fermi gas. 
The applicability of the mean-field theory in the weakly 
coupling or intermediate coupling regimes will be com- 
mented. In Sec. Ill, we will discuss in detail how the 
density profiles and the order parameters evolve with in- 
creasing temperatures. Together with the analysis of the 
entropy, energy and specific heat, we obtain a phase dia- 
gram at finite temperatures. Finally, Sec. IV is devoted 
to the conclusions and remarks. 
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II. SELF-CONSISTENT BOGOLIUBOV-DE 
GENNES THEORY 

Fermi gases of ^Li atoms near a broad Feshbach reso- 
nance can be well described using a single channel model, 
as confirmed both experimentally [gOl and theoretically 
(Hll, [i^l- To obtain the phase diagram in Fig. 1, we use 
the self-consistent BdG theory ^,i5l,[53|, by assuming a 
pairing order parameter A(a;) that breaks the U{1) sym- 
metry of the number conservation of total neutral atoms. 
This weak coupHng theory has been previously appHed to 
discuss the zero temperature properties of ID polarized 
Fermi gases by the present authors. We refer to Ref. [i^l 
for the detailed description of the model and the BdG 
formaHsm. We outHne below some essential ingredients 
of the theory. 

The BdG equations describing the quasiparticle wave 
functions u,, (x) and Vrj (r) , with excitation energies Ej) 
are [M]: 
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where the single particle Hamiltonian = 
-lT?V'^/2m -f gionwix) + Vextix), and Vext (x) = 
muj'^x'^/2 is the harmonic trap potential. The inter- 
action is a short range potential gmS^x), resulting in 
a diagonal Hatree term gionwix) and an off-diagonal 
order parameter potential in the mean-field decoupling. 
To account for the unequal spin population Na for the 
pseudospins a =t,J., the chemical potentials are shifted 
as = /i ± (5/i. This leads to different quasiparticle 
wave functions for the two spin components. However, 
there is a symmetry of the BdG equations under the 
replacement u*^{x) Vrj^{x), v^^^ix) - Urj^{x), 
Ejji ^^vP Thus, in Eq.(l) we have kept only the 
spin-up part, and obtained the solutions with both 
positive and negative excitation energies. 

The order parameter and the particle density of each 
component that appears in the BdG equations can be 
written as, 

A{x) = -giDY.u,{x)v;{x)f{E,), (2) 

n^ix) = ^<,(x)u,,(x) /(£;,,), (3) 
n 

'^i(^) = ^'"vix>vi^)fi-En), (4) 



where f{Ejj) = l/(exp[£'^/A;Br] -|- 1) is the Fermi distri- 
bution function. The particle density of each component 
is subjected to the constraint / dx (a;) = which 
eventually determines the chemical potentials. Eq. 
together with the definitions of the order parameter and 
particle densities, Eqs. ([2]), Q, and |[4|, form a closed 
set of equations, which has to be solved self-consistently. 
We have done so via a hybrid procedure, in which the 
high-lying excitation levels above an energy cutoff Ec has 



been solved approximately using a local density approxi- 
mation. This procedure is very efficient. We refer to Ref. 
[13, nil for further details. Note that the final numerical 
results are independent of the cutoff energy Ec, provided 
that it is large enough. 

Once the BdG equations are solved, it is straightfor- 
ward to calculate the total entropy and total energy of 
the system. The expression for the entropy is given by, 

S^-ksY. [/(^'O ^^-^fi^v) + fi-E„) In f{-E,)] , 

(5) 

where the summation can be restricted to energy lev- 
els with energy < E^. Other high- lying branches 
have an exponentially small contribution because of the 
large value of Ec{^ kBT). This has been neglected. 
In contrast, the total energy includes two parts, E = 
Edisc + Econt ■ The discrete contribution Edisc takes the 
form, 



nEr,)- / dx\v^{x)Y 



^l^N^+^ilNl~glD / dx\A{x)\ 

|B„l<Bc 

while the high-lying levels contribute 

Econt ^ J '^^ [^1 (^) + (x)] 



(6) 
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where ei (x) and 62 (x) are given by. 
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(9) 



respectively. An important self-consistency check of the 
calculations of the entropy and of the energy is given by 
the thermodynamic relation, 

E{T) = Eo + £dT'r(^^y (10) 

where Eq is the ground state energy at zero temperature. 

Before closing the section, we discuss briefiy the ap- 
plicablity of mean-field Bogoliubov-de Gennes theory in 
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one dimension. It is well-known that with decreasing di- 
mensionality, the pair fluctuation becomes increasingly 
important. Specifically to ID, the true long-range order 
is completely destroyed by fluctuations. Thus, strictly 
speaking, there is no a priori justiflcation for the use 
of mean-field BdG theory for a uniform spin-polarized 
Fermi gases. Nevertheless, we found in the previous 
study that at zero temperature and at a coupling con- 
stant 7 ~ 1, the mean-field BdG calculation provides 
reasonable description for the energy and chemical poten- 
tial of a uniform spin-polarized Fermi gas, as compared 
to the exact Bethe Ansatz solutions (See, for example, 
the Figs. 12 and 13 in Ref. 

On the other hand, the presence of a harmonic trap ef- 
fectively increases the dimensionality of the system. For 
instance, the density of state of the trapped ID system 
is a constant, exhibiting the same behavior as a uni- 
form Fermi gas in two dimensions. Thus, in a harmonic 
trap, we antipicate that the pair fiuctuations should be 
much suppressed, and thus the use of the mean-field BdG 
theory may be justified. We have checked in the previ- 
ous work the vaHdity of BdG theory for a trapped spin- 
polarized Fermi gas and found that indeed, at 7 ~ 1 the 
density profile of a trapped gas obtained from the zero 
temperature BdG calculations agrees well with that from 
the asymptotically exact Gaudin solutions (See, for ex- 
ample. Fig. 18 in Ref. [42]). 
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III. RESULTS AND DISCUSSIONS 

For concreteness, we consider a trapped polarized gas 
with a number of total atoms N = N-^ + Ni = 100. 
In traps, it is convenient to take trap units with, m = 
H — Lu — ks = 1, so that the length and energy will 
be measu red in units of the harmonic oscillator length 
aho — yh/muj and of the level spacing hw, respec- 
tively. The characteristic energy scale may be given by 
the Fermi energy of an unpolarized ideal gas at zero tem- 
perature, Ep = Nfiiojl^ which provides also a charac- 
teristic temperature scale, Tp = Ep/ks- On the other 
hand, the length scale is set by a Thomas-Fermi radius, 
xtf = N^/^aho- To characterize the interaction, we use 
a dimensionless coupHng constant at the trap center [47| . 
70 — aho/iT^N^/^aio), where aio — —mgiD/{nh^) is the 
ID scattering length, and 70 is roughly the ratio of the in- 
teraction energy density at the trap center to the kinetic 
energy density. Thus, 70 -C 1 corresponds to the weakly 
interacting Hmit, while the strong coupling regime is re- 
alized when 7o ^ 1. Throughout the paper, we use a 
coupling constant 70 = 1.6, c omp arable to the estimates 
for a real experimental setup [43|. 



A. Density profiles and order parameter 

Figs. 2 and 3 give BdG results for the density pro- 
files of each component, as well as the density difference 



Figure 2: (color online) Spin up and spin down density profiles 
(thick solid lines), the density difference (thin solid lines), and 
the order parameter (dot-dashed lines) of a trapped ID Fermi 
gas for a spin polarization P — 0.075, at several temperatures, 
T = O.OITf (upper panel), T = O.OSTf (middle panel), and 
T = O.ISTf (bottom panel). For a better illustration, the 
scale of the density difference has been doubled. At zero tem- 
perature, the system stays in a phase separation phase with 
a FFLO state at center and an outside BCS shell. The FFLO 
phase at the trap center, signaled by the oscillations in the or- 
der parameter, suppresses with increasing temperature, and 
leaves out a pure BCS state throughout the trap at about 
O.IOTf. At a higher temperature around 0.20rF, the system 
becomes fully normal through a second order superfluid phase 
transition. 



5n (x) — n| {x)—ni {x) and the order parameter A(a;), at 
different temperatures as indicated and at two spin po- 
larizations, P = {N^-N{)/N ^ 0.075 (Fig. 2) and 0.175 
(Fig. 3). In the work by Orso [45] and our previous 
studies [i^, 1131, it has shown that at zero temperature 
there are two phase separation states: For a small spin 
polarization the system stays in the spatially inhomoge- 
neous FFLO superfluid state at the trap center and in a 
BCS superfluid state at the edge (referred to as FFLO- 
BCS later). For a large spin polarization, the gas still 
remains in the FFLO state at the trap center, but be- 
comes a fully polarized normal state towards the edge of 
the trap (referred to as FFLO^N later). We refer to the 
Sees. VI and VIIA of Ref. [4'^ for further details. The 
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Figure 3: (color online) Same plots as in Fig. 2, but for 
a larger spin polarization P = 0.175. With this value of 
spin polarizations, at zero temperature the atoms at the trap 
edge stay in a normal state. With increasing temperature, 
the FFLO phase at the trap center vanishes gradually, which 
eventually leads to a normal phase throughout the trap at 
T ~ 0.15Tf. 




0.00 



0.01 



1/7V 



0.02 



Figure 4: (color online) (a) Density difference of a spin- 
polarized Fermi gas at P = 0.175 and T = O.OlTp for dif- 
ferent numbers of total atoms as indicated, (b) The number 
dependence of the amplitude of the density oscillations A at 
the trap center. The amplitude A decreases slowly with in- 
creasing the number of atoms, (c) The number dependence of 
the periodicity of the density oscillations d at the trap center. 



two spin polarizations in Figs. 2 and 3 are selected in 
such a way that initially (at low temperature) the cloud 
is in diflFerent phase separation ground state. We then 
trace how these two phase separation states develop as 
the temperature increases. 

For a small spin polarization below a critical value Pc, 
the FFLO-BCS phase separation state at low tempera- 
ture is clearly characterized by the spatial profile of the 
order parameter (Fig. 2a). It shows oscillations at the 
trap center, characteristic of FFLO states, and two shoul- 
ders at the trap edge, characteristic of a ID BCS state. 
Naively, the FFLO state is much more easily disturbed 
by thermal fluctuations than the BCS state, due to a re- 
duced pairing phase space. Therefore, as shown in Figs. 
2b and 2c, with increasing temperature the FFLO super- 
fluid is suppressed faster than the BCS shoulders at the 
edge. This leads to a pure BCS state above a certain 
temperature about O.lOTp. Further increase of the tem- 
perature will destroy the superfluid state eventually at 
around 0.20Tf. 

Experimentally, the change from FFLO-BCS phase 



separation state to the pure BCS state and flnally to the 
normal state may be monitored by measuring the density 
difference. At low temperature in FFLO-BCS separation 
phase, the local spin polarization or density difference is 
fully carried by the FFLO state and thus is restricted to 
the trap center (Fig. 2a). Here, we have neglected the 
single peak of density difference exactly at the trap edge 
that is caused by breakdown of the mean-fleld approach, 
see Ref. |47,] for detailed discussions. At a higher tem- 
perature, the BCS state starts to contribute to the spin 
polarization, due to the thermal excitations. As a result, 
the density difference leaks out gradually to the trap edge 
(Fig. 2b). When the FFLO state fully disappears, the 
difference in density becomes very flat throughout the 
trap (Fig. 2c). 

The temperature evolution for a FFLO-N phase sep- 
aration state with large spin polarizations is simpler, as 
shown in Fig. 3. In this case, the FFLO state at the 
trap center is destroyed gradually by increasing temper- 
ature, and one ends up with a fully normal cloud around 



6 




/» = 0.075 




A^=100 




0.58 



-1.0 -0.5 0.0 0.5 

X (N'^'a ) 



1.0 



Figure 5: (color online) Order parameters of a spin-polarized 
Fermi gas at P = 0.175 and T = O.OlTp for different numbers 
of total atoms: (a) = 50, (b) N = 100, and (c) iV = 200. 
The two dashlines indicate the value of the order parameter 
at the trap center. 



— ^=0 

- -P= 0.075 
P= 0.175 




Figure 6: (color online) Temperature dependence of the chem- 
ical potential (panel a) and of the chemical potential differ- 
ence (panel b) at three spin polarizations: P = (solid lines), 
P = 0.075 (dashed lines), and P = 0.175 (dash-dotted lines). 



of-mass momentum ^lo proportional to the Fermi wave- 
vector diflference, ^lo — kpi — kpi — {N^^^P)aJ^^ . Thus, 
in units of (N^^'^aho) the periodicity of the density oscil- 
lation (Fig. 4c) and of the order parameter are inversely 
proportional to the number of total atoms. 



T ~ O.lSTp. The shape of density difference profiles is 
nearly independent of temperature. Thus, the distinct 
temperature dependence of the density difference in the 
FFLO-BCS and FFLO-N states provides a useful way to 
distinguish these two phase separation phases. 

At the end of this subsection, we emphasize that in 
Figs. (2a) and (2b) it is the oscillation of the order pa- 
rameter that is the unambiguous signature of the FFLO 
state. There is also a related oscillation in the density 
profiles with doubling periodicity. However, these density 
oscillations are presumably due to finite size effect. We 
have checked this point by varying the number of total 
atoms. In Figs. 4 and 5, we show the density difference 
and the order parameter at a spin polarization P = 0.175 
and at a low temperature T = O.OlTp for particle num- 
bers in the range 50 ~ 200. With increasing the number 
of total atoms, the amplitude of the density oscillation at 
the trap center becomes smaller. In contrast, the ampli- 
tude of the order parameter A (a; = 0) stays nearly con- 
stant. Note that in Fig. 5 we find a Larkin-Ovchinnikov 
(LO) state in the center region, i.e., the order parameter 
behaves like A (x) = A (a; = 0) cos{qi^ox) with a center- 



B. Equation of state 

Figs. 6 and 7 show the equation of state as a func- 
tion of temperature. The total entropy and total energy 
increase monotonically with increasing temperature and 
spin polarization. In particular, at a finite spin polariza- 
tion, the entropy increases linearly with temperature, as 
compared to the exponential temperature dependence for 
a BCS superfiuid {i.e., P — 0). This linear dependence 
arises from having gapless single-particle excitations lo- 
cated at the node of FFLO states at the trap center, and 
in case of the FFLO-N phase, from the normal cloud at 
the edge of the trap. As the temperature increases, non- 
trivial kinks are visible in both chemical potentials and 
entropy, refiecting the smooth transitions between differ- 
ent states as mentioned earlier. 



C. Specific heat and finite temperature phase 
diagram 

The analysis of the profiles of the density distributions 
and of the order parameter has already given some qual- 
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Figure 7: (color online) Temperature dependence of the total 
entropy per particle (panel a) and of the total energy per 
particle (panel b) at three spin polarizations: P = (solid 
lines), P = 0.075 (dashed lines), and P = 0.175 (dash-dotted 
lines). 
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Figure 8: (color online) Temperature dependence of the spe- 
cific heat at three spin polarizations: P = (solid lines), 
P = 0.075 (dashed lines), and P = 0.175 (dash-dotted lines). 
For P — 0.075, the dashed arrow indicates the transition from 
a FFLO-BCS phase separation phase to a pure BCS state, 
while the solid arrow marks the second order phase transition 
from a BCS superfiuid to a normal gas. 



itative features for the phase diagram at finite tempera- 
tures. To quantify it, we calculate the specific heat, using 
the definition, 

_dE_ dS 



Fig. 8 displays the temperature dependence of the spe- 
cific heat at three spin polarizations. Different transitions 
are well characterized by the apparent kinks in the spe- 
cific heat. For a small polarization, two kinks are clearly 
visible, corresponding to the transition from a FFLO- 
BCS phase separation state to a pure BCS state, and 
in turn to a normal state. In contrast, for a large po- 
larization, typically only one kink is identifiable, which 
should be attributed to the transition from a FFLO-N 
phase separation state to a completely normal state. 

Gathering the position (temperature) of the kinks 
for different spin polarizations, we obtain the solid and 
dashed lines in the phase diagram, as shown in Fig. 1. 
The former line corresponds to the second-order phase 
transition from either a BCS or FFLO-N to the normal 
state; while the later distinguishes a high-T pure BCS 
state from a low-T FFLO-BCS phase separation state. 
The determination of the boundary between FFLO-BCS 
and FFLO-N phase separation phases, for example, the 
critical spin polarization Pc{T), is more difficult from 
numerics. As discussed in our previous works [i^ . [47| . 
the critical spin polarization is calculated from a criti- 
cal chemical potential difference, i.e., the half of binding 
energy of the ID molecule pairs. From Fig. 6b, for a 
large spin polarization the chemical potential difference 
is temperature insensitive. On quite general grounds, we 
thus assume that the critical spin polarization Pc{T) is 
nearly temperature independent. This gives the dash- 
dotted line in the phase diagram (Fig. 1). 

It is evident from Fig. 1 that there is a wide tempera- 
ture window for the presence of FFLO states at the trap 
center. The typical temperature for observing the FFLO 
state would be of one-tenth of the Fermi temperature, 
corresponding to an entropy S ^ 0.2NkB. Such a tem- 
perature or entropy is within the reach of present-day 
techniques 



IV. CONCLUDING REMARKS 

In conclusion, based on a mean-field Bogoliubov-de 
Gennes theory, we calculate the thermodynamic proper- 
ties of a polarized atomic Fermi gas in a highly elongated 
harmonic trap. The profiles of the density distributions 
and of order parameter, the equation of state, as well as 
the specific heat have been analyzed in detail. We have 
then established a finite temperature phase diagram. Our 
results are useful for experiments at Rice University [H3|, 
in which a search for the exotic FFLO states in ID po- 
larized Fermi gases is under-taken. Our estimated tem- 
perature and entropy for the realization of FFLO states 
are about O.lTp and 0.2iVfc_B, respectively. These values 
are experimentally attainable [7]. 

While our weak-coupling Bogoliubov-de Gennes study 
provides a semi-quantitative finite-temperature phase di- 
agram of the ID polarized Fermi gas, an improved de- 
scription could be obtained by solving the exact thermo- 
dynamic Bethe ansatz solutions of the ID gas [H^l, with 
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the trap eflFect treated using a local density approxima- 
tion [43, [4fl|. We plan on doing this in future work. 
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